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1. Introduction 

The concept of entanglement is intimately related to the mathematical structure of 
quantum mechanics, particularly as a direct consequence of linearity in tensor product 
Hilbert spaces. Apart from its theoretical relevance as one (if not the most) striking 
features of quantum mechanics, entanglement has been shown to enable (and play a 
crucial role in) practical applications such as teleportation protocols and superdense 
coding. It is, for this reason, often regarded as a resource which needs to be studied 
from several standpoints. 

A considerable amount of research has been devoted to unveil the mathematical 
structures underlying entanglement, in particular concerning its quantification. 
Properties for a good entanglement measure are reviewed and discussed in [8j [9j [15] . An 
alternative to the analytic approach was presented by Brown et al [I] and later by Borras 
et al [lj. Both works explore the application of numerical optimisation techniques to 
search for highly entangled multi-qubit states. 

In the work of Brown et al [I], the search is applied over the set of mixed states 
using the negativity of the partial transpose as a measure of entanglement. The search 
procedure is basically a hill-climbing algorithm with some adjustments. Such a simple 
form of search suffices, as the cost function which is sought to maximise (i.e. the 
entanglement as defined by the negativity) is known to be convex. Among the results 
reported are the following two 5-qubit states 

1*")= 2^(|00110) + |01011) + |10001) + |11100) 

+i(|00101) + |01000) + |10010) + 1 11111))) ^ ' 

1*"')= 2^(|00110) + |01001) + |10101) + |11010) 

+i(|00000) + |01111) + |10011) + |11100))) ^ ' 

which are essentially equivalent and can be reformulated as 

= ^(|001>|$_) + |010)|¥_) + 1100)1$+) + |lll)|¥ + » (3) 

where |*±) = ^=|00) ± |11) and |$±) = ^=|01) ± |10). The state \i/> 5 ) achieves the 
maximum possible entanglement according to the negativity measure. Muralidharan 
and Panigrahi (T2J [13] have recently investigated the usefulness of this state for quantum 
teleportation, superdense coding and quantum secret sharing. In particular, perfect 
teleportation of arbitrary 1- and 2-qubit systems is possible, as well as superdense coding 
of 5 classical bits by using three qubits. (Incidentally, it is also shown how to physically 
produce this state in a relatively easy manner.) 

In the work of Borras et al [lj, the search is restricted to pure states of 4, 5, 6 
and 7 qubits, and several measures of entanglement are investigated. The overall search 
algorithm is very similar to the hill-climbing scheme employed by Brown et al , and the 
authors also explore the statistical distribution of entanglement for systems of the sizes 
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aforementioned. In the case of 6 qubits, they discovered the state 



iv*> = w 



loooooo) + 

|111010) 
1001111) + 
|101101) + 
1001010) 
|101011) + 
|011110) + 



1000011) + 
|111001) + 
1010001) + 
|100111) + 
1001100) + 
1 101000) + 



|111100) + 
1001001) + 
|101110) + 
|011101) + 
|110011) + 
|011011) + 



|000101) + 
|110110) + 
|010010) + 
1 100010)- (4) 
|010100) + 
|100100) + 



1111111) + 

1000110) + 
1 110000) + 
1011000) + 
|110101) + 
|010111) + 
jlOOOOl)) 

which, as in the case of Brown et al 's \ips), possesses a rather simple structure and 
maximum entanglement according to several measures. Choudhury et al [5] have also 
found various applications for this state, namely quantum teleportation of an arbitrary 
3-qubit state and quantum state sharing of an arbitrary 2-qubit state. 

In the case of 7 qubits, the authors report a highly entangled state (approximately 
155.812856, according to the negativity measure) with, apparently, no simple structure. 
The state is given as a list of 128 complex numbers corresponding to the coefficients of 
the superposition, and has all the 1-qubit density matrices completely mixed, but not 
the 2- and 3-qubit ones. 



1.1. Motivation 

Apart from theoretical repercussion, genuinely entangled states of 5 and 6 qubits have 
found practical applications and new protocols based on them have been derived. In the 
case of 7 qubits, it is not known how to write "nicely" the state (s) reported by Borras 
et al . This lack of structure may be a drawback in deriving protocols based on it. 

The search schemes described above have proven to be very successful, but the 
resulting states can be difficult to process manually. In this work, we present an 
alternative formulation of the search where not only high entanglement is sought, but 
also some desired structure is imposed on the form of the outcome. As a proof-of- 
concept, we show how highly entangled states with a very compact form can be easily 
found. 



2. Computationally feasible entanglement measures 

All the numerical optimisation approaches aforementioned rely directly upon the use of 
a computationally tractable measure of entanglement. For completeness and readability, 
we next provide a brief introduction to some of them, and particularly to the measure 
used in this paper. For a more in-depth treatment of the subject, we refer the reader to 

The degree of bipartite entanglement of a composite quantum system can be 
quantified in terms of the purity (or mixedness) of the reduced density matrix of one 
of the two subsystems: the lower the purity (or the higher the mixedness), the higher 
the entanglement. There is not, however, a unique way of extending this notion to 
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systems composed of more than two subsystems. A natural alternative is to quantify 
multipartite entanglement among n parties as a function of the mixedness of all the 
possible bipartitions of the whole system. A common, though not unique, approach is 
to take the sum of the individual entanglement measures associated with the 2"~ 1 — 1 
possible bipartitions of a n-qubit system [U [1] 



where p is the density matrix of the n-qubit system, p s is the reduced density matrix 
associated with subsystem s, and the sum is taken over all the possible bipartions. 
Common examples of bipartite entanglement measures Eb are: 

(i) The linear entropy: Sl(p) = 1 — Tr(p 2 ) 

(ii) The von Neumann entropy: Svn(p) = — Tr(plog 2 p) 

(iii) The Renyi entropy with q — > oo: S^°°(p) = — lnA^ ax , Xk being the eigenvalues of 



All these entropic measures constitute clear examples of computationally feasible 
measures of entanglement for pure states. They cannot, however, distinguish between 
classical and quantum correlations, and so their use for mixed states is quite limited. 

Brown et al developed in [I] an entanglement measure for mixed states based on the 
negative partial transpose (NPT) criterion, which establishes as a necessary condition 
for separability of any density matrix that it has only non-negative eigenvalues [14J. 
Details about how to compute the partial transpose can be found in [3]. 

The negativity, first introduced in [17] and subsequently explored in [161 Ell IH EE], 
is defined as the sum of all the negative eigenvalues of the 2 n ~ 1 — 1 partial transposes. 
The result if often negated in order to have a positive measure. We will refer to this 
measure as E^p^. In some cases it is convenient to have a measure normalised between 
and 1, which can be achieved by dividing by the maximum possible entanglement. 
Upper bounds for the four entanglement measures discussed above can be derived by 
considering a hypothetical n-qubit pure state such that all its marginal density matrices 
are completely mixed pQ. In the case of the £"npt 5 this maximum is 1.5, 6.5, 17.5, 60.5, 
157.5, 504.5, 1297.5, ... for systems composed of 3, 4, 5, 6, 7, 8, 9, . . . qubits. We will 
denote by i^pr the normalised entanglement. 

3. Searching for simplicity 

3.1. Search space 

We want to search for pure states having high entanglement and simple algebraic 
structure. There is not, however, a unique way to define what "simple" means, and 
some of such definitions might not be computationally appropriate. For the purposes of 




(5) 



s 
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Table 1. Upper bounds for |§"(V)| for some values of n and |V|. 



Number of qubits 



\v\ 


4 


5 


6 


7 


8 


5 


2 37.1 


2 74.3 


2148.6 


2297.2 


2594.4 


7 


244.9 


289.8 


2 179.7 


2359.3 


2 718.7 


9 


2 50.7 


2101.4 


2202.9 


2405.7 


2811.5 


11 


255.3 


2 110.7 


2221.4 


2442.8 


2885.6 


13 


2 59.2 


2 118.4 


2236.8 


2473.6 


2947.3 



this work, simple states will be those which can be written in a "nice" way with respect 
to a given basis, i.e. states of the form 

2 n -l 

IV>> = ^ £ *l*> (6) 

i=0 

where only a few of the coefficients q are non-null and, in turn, are nice to write 
(e.g. ±1, ±i, ±(1 ± i), . . .). Even though this is quite an arbitrary definition of 
simplicity, it captures well the intuitive idea of what a nice-to-write state means; and, 
more importantly, provides a measurable way to determine how simple a state is, a 
crucial property if we are to somehow incorporate such an aspect into the search. 
For the purposes of this paper, the space of "simple" states will be given by 

2"-l 

where 1/K is the appropriate normalisation factor. The size and properties of such a 
spacdj] strongly depend on the set V of allowed coefficients. Three examples that will 
be used later are V 3 = {0, ±1}, V 5 = {0, ±1, ±i} and V 9 = {0, ±1, ±i, ±(1 ± i)}. 

Using a simple counting argument it is not difficult to see that the size of the search 
space is |§ n (V)| = \V\ 2n . Note, however, that some states might actually be identical 
up to some global phase, and therefore previous expression is actually an upper bound 
on the number of useful states. The exact amount of indistinguishable states depends 
on the particular elements in the set V and we will not be generally concerned about 
them. Table Q] shows the size of the search space for the number of qubits considered 
in this work and for some small values of \V\. The combinatorial explosion is easily 
recognisable in these figures, and any attempt of search by enumeration is beyond our 
current computational resources even for small values. 

3.2. Fitness function 

A good fitness (or cost) function should provide guidance through the search space 
towards the desired solutions, regardless of whether or not its numerical values are 

\ Note that "space" is used here as in "search space" and not as in "Hubert space", as S n (V) might 
not always be one. 
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related to the actual cost or benefit of particular solutions. As we are interested in 
states having high entanglement and as many null coefficients as possible, a natural 
choice for the fitness function is to explicitly reward both features. We have used a 
fitness function of the form 



and a G [0,1] is a punishing factor. The term E^p™\ip) (ip\ simply computes the 
entanglement as measured by the NPT, while the right part counts the percentage 
of non-null coefficients in the superposition. Note that both terms are normalised in the 
interval [0, 1], so a admits a natural interpretation as how much punishment is placed 
on the amount on non-null coefficients. Other variations around the same idea are of 
course possible. 

3. 3. Move function 

The move function determines which states are reachable from a given one at any point 
of the search. We have defined a move function consisting of randomly choosing one 
of the 2" coefficients of the state, q, and another, Cj, from V such that c, ^ Cj. The 
coefficient q is replaced by Cj and the state vector is then renormalised. Note that each 
state vector has exactly 2™(|l / | — 1) different^ neighbours. We will write \<p) <— Move|?/>) 
to denote that state \<p) is the result of moving state \ip). 

3.4- Search procedure 

It is well known that the negativity cost function is convex [16] and its smoothness 
provides us with a relatively easy-to-walk search landscape. In particular, its convexity 
for mixed states implies that there is no chance for the search of being caught in a 
local optimum, even though the case for pure states is less clear. However, the trade-off 
between entanglement and simplicity contained in expression ([8]) translates into a more 
rugged search landscape, with lower correlation between the fitness of neighbour state 
vectors and, therefore, a substantial increase in the number and distribution of local 
optima. This suggests that very simple forms of search as those used in previous works 
(basically hill-climbing algorithms) are not likely to be appropriate in this case. We have 
verified experimentally this fact and concluded that a more powerful search technique 
is necessary. 

Simulated Annealing [11] is a well-known search heuristic inspired by the cooling 
processes of molten metals. It can be seen as a basic hill-climbing coupled with the 

§ As discussed above, some of them might be identical up to a global phase. 
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probabilistic acceptance of non-improving solutions. This mechanism allows a local 
search that eventually can escape from local optima. 

The search starts at some initial state (solution) \ip) Er S n . The algorithm employs 
a control parameter T G R + known as the "temperature" . This starts at some positive 
value T and is gradually lowered at each iteration, typically by geometric cooling: 
Ti+i = (3Ti, p G [0, 1]. At each temperature, a number MIL (Moves in Inner Loop) of 
neighbour states are attempted. A candidate state \<f>) in the neighbourhood of \ip) is 
obtained by applying the move function. The new state is always accepted if it is better 
than To escape from local optima, the technique also accepts candidates which are 
slightly worse than the current state, meaning that its fitness is no more than |Tln[/| 
lower, with U a uniform random variable in (0, 1). As T is gradually lowered, this term 
gets closer to and it becomes harder to accept worse moves. 

The algorithm terminates when some stopping criterion is met, usually after ex- 
ecuting a fixed number of inner loops or when some maximum number of consecutive 
inner loops without improvements have been reached. The basic algorithm is shown 
below. 

1 m e R S n 

2 \*p max ) <- 

3 T <- T 

4 repeat until stopping criterion is met 



5 repeat MIL times 

6 \</>) <- Move|-0) 

7 U £ R (0, 1) 

8 if F\<p) > F\ip) +T\nU then 

9 «- |0) 

10 if F\4>) > F\4> max ) then 

11 |^ma^ ^_ |^ 

12 endif 

13 endif 

14 endrepeat 

15 T <— (5T 



16 endrepeat 

17 return \if; max ) 

3.5. Implementation 

The search scheme described above was implemented in C++, using the GNU 
Scientific Library [6] for complex numbers, matrix manipulations, pseudorandom 
number generation, eigenvalues, etc. 

4. Results 

We next present and discuss some results obtained for systems of 4, 5, 6, 7 and 8 qubits. 
A couple of states are shown for each case. These are representative of the typical 
outcome of the search, but many others have been produced by the technique. 
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In all the cases, the search algorithm is parameterised with 1000 moves in the 
internal loop (MIL) and a stop criterion consisting of 10 consecutive internal loops 
without improvement. For the number of qubits studied, the best results are achieved 
with a between 0.1 and 0.3, a starting temperature T between 0.0005 and 0.001, and a 
cooling rate (3 between 0.99 and 0.999. 



4-1. 4-qubit states 

Higuchi and Sudbery proved in [7j that there is no 4-qubit pure state having all its 
marginal density matrices maximally mixed. Consequently, the theoretically maximum 
amount of entanglement for systems of 4 qubits (e.g., 6.5 according to the negativity 
measure) is unachievable. The state 

\HS) = ^(|1100) + |0011)+cj(|1001) + |0110))+cj 2 (|1010) + |0101)))(10) 
V6 

where u = — ~ + ^i is also reported in [7j. It is conjectured to be maximally entangled, 
having Enpt\HS) ~ 6.0981. Borras et al claim in [1] that their search procedure 
consistently finds states with exactly this same amount of entanglement. In fact, the 
HS state is known to be a local maximum for the von Neumann entropy [3], and some 
arguments given in [7] suggest that it may also been a global maximum for this measure. 
In our experimentation using the set of coefficients V5 we found the state 

K) = 4=(|1001) - z(|0000) - |0011) + |0110) + |1100) + |1111))) (11) 
V6 

EnptK) « 5.989631 (12) 

This state has 22 negative eigenvalues from 7 partial transposes. For the four single- 
index cuts, each partial transpose has a single negative eigenvalue of —5, and all the 

1- qubit density matrices are maximally mixed. The 18 two-index partial transposes 
have a more complex structure, each one having 6 negative eigenvalues. None of the 

2- qubit density matrices are maximally mixed. 

When using the set V9 we found the state 

M = ^(10001) - 10100) H-»(|1011> - |1110))+ 

(1 + z)(|0010) + 1 1101) - |0101) - |0110) - |1010) - |1100))+ (13) 
(1 - i)(|1000) + |1001) - 10011) - |0111))) 

E NPT \i)l) « 6.051660 (14) 

As in the case of \tp2), this state has 22 negative eigenvalues from 7 partial transposes. 
The four single-index cuts give a partial transpose with a single negative eigenvalue 
of —1, and all the 1-qubit density matrices are maximally mixed. The 18 two-index 
partial transposes have 6 negative eigenvalues each one. Again, none of the 2-qubit 
density matrices are maximally mixed. 

The HS state is obviously more entangled than \ifjf) and l^), yet both states 
(particularly the latter) do have a very high entanglement and a rather simple structure. 
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Merely for validation purposes, we also attempted searches including both u and uj 2 in 
the set of valid coefficients V. In such cases, the HS state (or states essentially equivalent 
to it) is easily found. 

4-2. 5-qubit states 

The search successfully finds 5-qubit states exhibiting maximum entanglement (-Enpt = 
17.5) and a rather simple algebraic structure. The following two, and are 
particular examples found using coefficients in V5 and Vg, respectively 

|Vf> = 2^(101010) + 110011) + 110110) + 111000) 

+z(|00001) - 100100) + 101111) - 1 11101))) [ ' 

£npt#5> = 17 - 5 ( 16 ) 

\4) = + O(ionoi) + 101110) - 110100) + 110111) 

+ |11000) + |11011)) + (1-«)(|00001) - 100010))) 1 ' 

^nptIV's) = 17-5 (18) 

Both states have 65 negative eigenvalues from 15 partial transposes. For the five 
single- index cuts, each partial transpose has a single negative eigenvalue of — |. For the 
ten two- index cuts, each partial transpose has six negative eigenvalues of — \. So 

£ NPT |V>5> = £npt|^> = 5 • \ + 10 • 6 • 1 = 17.5 (19) 

All the partial density matrices are completely mixed for both states, i.e. Tr(p 2 ) = | 
for all the 1-qubit marginal density matrices, and Tr(p 2 ) = | for all the 2-qubit marginal 
density matrices. 

It is not difficult to show that these two states are equivalent under local unitary 
transformations, and also equivalent to Brown et al 's ^5) state given by ([3]). All the 
maximally entangled 5-qubit states found in our experimentation have this form. 

4-3. 6-qubit states 

As with the case of 5 qubits, the search finds maximally entangled (-Enpt = 60.5) 6-qubit 
states having simple structure. Two examples are provided below 

|$») = 1( - 1000101) - 1001011) - 1010001) + 1011111) 
+ 1100000) + 1101110) - 1110100) + 1111010) 

+z(|000110) - 1001000) - 1010010) - 1011100) 

+ 1100011) - 1101101) + 1 110111) + 1111001))) 
E NPT m = 60.5 (21) 

|^) = i(|011010) - 1011111) + 1101011) - 1101110) 
-1110011) - 1 110110) - 1111000) - 1111101) 
+i(-|000010) - 1000111) + 1001001) + |001100) 
+ 1010001) - 1010100) - liooooo) + 1100101))) 



(20) 



(22) 
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£npt#6> = 60 - 5 ( 23 ) 

Both states have 376 negative eigenvalues from 31 partial transposes. For the 6 
single-index cuts, each partial transpose has a single negative eigenvalue of — |. For the 
15 two- index cuts, each partial transpose has 6 negative eigenvalues of —j. For the 10 
three-index cuts, each partial transpose has 28 negative eigenvalues of — |. So 

E NPT |<) = £ NP i#6> = 6 • ~ + 15 • 6 • ~ + 10 • 28 • i = 60.5 (24) 

All the partial density matrices are completely mixed for both states, i.e. Tr(p 2 ) 
being |, | and | for all the 1-qubit, 2-qubit and 3-qubit marginal density matrices, 
respectively. 



4-3.1. A new 6-qubit state The states given by (12"U1) and (12"2"1) are equivalent under local 
unitary transformations and can be reformulated as 

|* 6 > = ^(|Fo}|*_) + |F!>|* + > + |F 2 )|$_) + |F 3 >|$ + » (25) 

with |F ) = ^(|0000) + |1111)), \F X ) = ^(|0011) + |1100)), \F 2 ) = ^(|0110) + |1001)), 
and |F 3 ) = ^(|0101) + 1 1010) ) . Note that this state has an elegant description in terms 
of Bell pairs, since 

^o> + i^>)H* + >l* + > (26) 

^(|F 2 ) + |F 3 )) = |$ + )|$ + ) 1 > 

This state has arguably a simpler algebraic structure than Borras et al 's 1^) 
given by (jlj). We postpone for future work a further study of the possible relationships 
between both states. 



4-4- 7-qubit states 

It has been pointed out that entanglement in 7-qubit systems might exhibit some 
similarities with the case of 4 qubits. A conjecture by Borras et al p] establishes that 
there is no pure state of 7 qubits whose marginal density matrices for subsystems of 1, 
2, or 3 qubits are all completely mixed. If such a state would exist, its entanglement 
as measured by the NPT would be 157.5. In [lj it is reported that 7-qubit states 
with entanglement up to -Enpt ~ 155.812856 are found. They all have completely 
mixed single-qubit marginal density matrices, but not completely mixed 2- and 3-qubit 
marginal density matrices. 

The results found with our approach are consistent with this numerical evidence. 
The states and l^) shown below have high (but not maximum) entanglement. 
Contrarily to the 7-qubit state reported in [T], these states have a remarkably simple 
structure and, as will be discussed later, totally mixed 1- and 2-qubit marginal density 
matrices. 
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The state 

\^) = ^=((1 + i)(|0000010) + |1000101) + 11011011) - 10011001) 

-|1010000) - 11111101)) + (1 - i)(|0001100) + 10010111) , s 
+|0100001) + 10101111) + 10110100) + 10111010) + 11001110) 
+ |1100011) - |1101000) - 11110110))) 

£ NPT |^) = 152.646039 (28) 

have 1113 negative eigenvalues from 63 partial transposes. For the 7 single-index cuts, 
each partial transpose has a single negative eigenvalue of — \. For the 21 two-index cuts, 
each partial transpose has 6 negative eigenvalues of — |. For 21 out of the 35 three-index 
cuts, each partial transpose has 28 negative eigenvalues of — |; the remainder 14 partial 
transposes have each 28 negative eigenvalues, but different from — |. 

Both the 7 single-qubit and the 21 2-qubit marginal density matrices are all 
completely mixed, i.e. with Tr(p 2 ) being | and |, respectively. For the case of the 
35 3-qubit marginal density matrices, 21 of them are completely mixed, corresponding 
to the following bipartitions 

{0, 1, 2}, {0, 1,4}, {0, 1, 6}, {0, 2, 3}, {0, 2, 4}, {0, 2, 5}, {0, 2, 6}, 

{0,3, 5}, {0,4, 6}, {1,2, 4}, {1,2,6}, {1,3, 4}, {1,3,5}, {1,3, 6}, (29) 

{1,4, 5}, {1,5, 6}, {2, 3, 5}, {2, 4,6}, {3,4, 6}, {3, 5,6}, {4, 5, 6} 

The 14 remaining marginal density matrices associated with bipartitions 

{0, 1, 3}, {0, 1, 5}, {0, 3, 4}, {0, 3, 6}, {0, 4, 5}, {0, 5, 6}, {1, 2, 3}, 

{1,2, 5}, {1,4, 6}, {2, 3, 4}, {2, 3,6}, {2,4, 5}, {2, 5,6}, {3,4, 5} 1 ] 

are not completely mixed, having Tr(p 2 ) = 0.156250 in all the cases except for {1, 4, 6}, 
which is 0.218750. 

A different 7-qubit state with high entanglement and very simple algebraic structure 

is 

\$) = ^((1 + i)(|0011101) + 10100010) + 10111000) + 11101100) 

+ |1111011) - 10001001) - 11000111) - 11110101)) ( s 

+(l-i)(|0010011) + |0110110) + |1001010) + |1010000) ^ ' 
+ |1011110) - 10000100) - 10101111) - 11100001))) 

£ NPT |^) = 152.504073 (32) 

The analysis for \^\) is completely similar to the previous case. Again, all the 
marginal density matrices for subsystems of 1 and 2 qubits are completely mixed. 
In the case of 3-qubit subsystems, the 18 marginal density matrices corresponding to 
bipartitions 

{0, 1, 2}, {0, 1, 3}, {0, 1, 4}, {0, 1, 5}, {0, 2, 6}, {0, 3, 5}, 

{0,4, 6}, {1,2, 6}, {1,3, 5}, {1,4,6}, {2, 3, 4}, {2, 3,5}, (33) 

{2, 3, 6}, {2, 4, 5}, {2, 5, 6}, {3, 4, 5}, {3, 4, 6}, {4, 5, 6} 
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are completely mixed, though not the remainder 17 corresponding to bipartitions 

{0, 1, 6}, {0, 2, 3}, {0, 2, 4}, {0, 2, 5}, {0, 3, 4}, {0, 3, 6}, 

{0,4, 5}, {0,5, 6}, {1,2, 3}, {1,2,4}, {1,2, 5}, {1,3,4}, (34) 

{1,3,6},{1,4,5},{1,5,6},{2,4,6},{3,5,6} 

The squared marginal density matrix have in all cases a trace of 0.156250, except for 

the case {2,4,6}, which is 0.187500. 

4-5. 8-qubit states 

The search for 8-qubit states is not as successful as in previous cases. An important 
factor in this case is definitely the computational hardness of exploring numerically 
the space. Each evaluation of the fitness function on a candidate state takes around 
7 seconds. This constitutes a serious limitation to the amount of states that can be 
explored in reasonable time. With the computing power available to us, a typical 4-day 
search visits around 40000 states (apart from the fitness evaluation, the search algorithm 
also imposes an overhead). In the case of 5, 6, and 7 qubits the search finishes well 
before 40000 states are visited, but this seems not to be the case here. We empirically 
found that a possible reason for this poor performance is related to the size of the 
neighbourhood. According to our move function, each candidate state has 2 n (|l / | — 1) 
neighbours. For a fixed V, the size of the neighbourhood increases exponentially in the 
number of qubits. In the case of 8 qubits, the average number of neighbour states visited 
before a new state is accepted is considerably higher than in previous cases, resulting in 
a very slow progression of the search. 

According to the NPT measure, the maximum possible (attainable or not) 
entanglement for a 8-qubit state is 504.5. We have found states with entanglement 
up to around 440 and rather simple structure. As an example, the state 

life>= §( 

100000010) - |00000101) - 100001111) + 100010111) + 100011010)- 
100011011) - (00100001) - |00100101) + 100100110) + |00101000) + 
100101110) + 100101111) - 100110000) - 100110011) + 100110100)- 
100111100) - (01000000) + 101000110) + 101000111) - |01001010) + 
101001100) - 101001101) + |01010001) + 101010110) + |01011010) + 
101100010) - 101100011) + 101101001) + 101111001) + |01111101) + 
101111110) + (10000001) + 110000111) + (10001000) + |10001010)+ 
110001100) - 110001110) + (10010000) + |10010100) + |10011101) + 
|10100010) - |10101001) + |10101011) - 110110111) + 110111001)- 
110111010) - 110111111) - 111000011) - |11000101) - |11001011)- 
|11010000) + 111010011) + |11010100) - 1 11011011) + |11011111) + 
111100000) + 111100100) - 111101101) - 111110101) - |11110111) + 
111111000) + 111111010) - 111111100) + 111111110)) 

has -EnptIV^) = 439.302328 and only 64 non-null coefficients equal to ±1. The 8 
single-qubit marginal density matrices are almost completely mixed, with Tr(p 2 ) being 
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0.500488 (5 times), 0.500977 (3 times) and 0.5 (1 times). Something similar happens for 
the 28 2-qubit density matrices, the average Tr(p 2 ) being 0.253540. The 3- and 4-qubit 
marginal density matrices deviates considerably more, on average, from the expected 
value of a completely mixed state. 

4-6. Performance and state evolution 

For 5 and 6 qubits, most runs result in a maximally entangled state whenever the 
parameterisation used for the search algorithm is appropriate. In both solution 
is often found in the first 4-6 cycles (i.e. the search ends after 4000-6000 generations). 
This constitutes a fairly rapid convergence, a typical search taking around 20-30 seconds 
for 5 qubits and 1-2 minutes for 6 qubits on a usual laptop. The cases of 4, 7 and 8 
qubits are slightly different, as the search never converges to a maximally entangled state 
and the procedure ends after a certain number of consecutive non-improving cycles. A 
typical search resulting in a 7-qubit state with -Enpt > 152 takes around 3-4 hours, while 
the 8-qubit case requires around 4-5 days. The running time is obviously exponential 
in the number of qubits, and the fundamental bottleneck comes from the evaluation of 
the fitness function, particularly the entanglement measure. 

Figure [T] shows the evolution of a typical search for 7 qubits. Apart from the 
fitness, both entanglement and the fraction of non-null coefficients are shown. It can 
be clearly observed the efect of the punishment factor in the fitness function: after a 
few thousands evaluations, candidate states have approximately 80% of null coefficients. 
The overall behaviour is identical for searches of 4, 5, 6, and 8 qubits. (In the case of 5 
and 6, however, the entanglement component eventually reaches the upper bound and 
the search stops.) 

5. Conclusions and future work 

In this work we have shown how states with simple structure can be discovered by 
restricting the search to certain subspaces of interest. For 5 qubits, states equivalent 
to Brown et al 's ^5) are easily found. In the case of 6 qubits, a new state with an 
algebraic structure arguably simpler than Borras et al 's [ipo) state has been discovered. 
Similarly to what has been pointed out for the latter, new teleportation and secret state 
sharing applications could be constructed using them. Furthermore, the fact of having 
at our disposal many of such states might enable us to devise new applications based 
upon the use of multiple, non-identical states having maximum entanglement. 

For the case of 7 qubits, various states with high entanglement and very simple 
algebraic structure have been found. Even though the entanglement is not maximum, 
all the single- and 2-qubit marginal density matrices are completely mixed, as well as 
more than half of the 3-qubit marginal density matrices. These features make of them 
potentially useful to support some applications. 

The case of 8 qubits is hard to explore numerically with our available computational 
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Figure 1. Typical state evolution (in colour in the electronic version) 

resources. Although the experimentation has not been as intensive as before (and 
therefore the parameterisation employed might not be adequate for this case), our 
preliminary results seems promising. Compact 8-qubit states with higher entanglement 
might be found with more computational resources. 

In this work we have restricted the search to pure multi-qubit states. Consequently, 
the use of the negativity of the partial transpose as measure of entanglement is a valid 
but unnecessary choice. An entropic measure such as the von Neumann entropy should 
be enough in terms of providing adequate search guidance. 

An interesting avenue for future research is related to a further work of Borras 
et al [2j. Here the authors explore the robustness of entangled states under various 
decoherence channels. Along these lines, we plan to study the entanglement decay of 
the states introduced in this paper. 
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